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This paper briefly summarizes previous work on complex classical mechanics and 
its relation to quantum mechanics. It then introduces a previously unstudied area 
of research involving the complex particle trajectories associated with elhptic poten- 
tials. 
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I. INTRODUCTION 

In generalizing the real number system to the complex number system one loses the 
ordering property. (The inequality Zi < Z2 is meaningless if zi, Z2 G C.) However, extending 
real analysis into the complex domain is extremely useful because it makes easily accessible 
many of the subtle features and concepts of real mathematics. For example, complex analysis 
can be used to derive the fundamental theorem of algebra in just a few lines. (The proof 
of this theorem using real analysis alone is long and difficult.) Complex analysis explains 
why the Taylor series for the function f{x) = 1/(1 + a;^) converges in the finite domain 
— 1 < 2; < 1 even though f{x) is smooth for all real x. 

This paper describes and summarizes our ongoing research program to extend conven- 
tional classical mechanics to complex classical mechanics. In our exploration of the nature of 
complex classical mechanics we have examined a large class of analytic potentials and have 
discovered some remarkable phenomena, namely, that these systems can exhibit behavior 
that one would normally expect to be displayed only by quantum-mechanical systems. In 
particular, in our numerical studies we have found that complex classical systems can exhibit 
tunneling-like behavior. 

Among the complex potentials we have studied are periodic potentials, and we have 
discovered the surprising result that such classical potentials can have band structure. Our 
work on periodic potentials naturally leads us to study the complex classical mechanics 
of doubly-periodic potentials. We have chosen to examine elliptic potentials because such 
functions are analytic and thus have well-defined complex continuations. This is a rich and 
previously unexplored area of mathematical physics, and in the current paper we report 
some new discoveries. 

This paper is organized as follows: In Sec. |TT]we give a brief review of complex classical 
mechanics focusing on the complex behavior of a classical particle in a periodic potential. 



In Sec. III we describe the motion of a particle in an elliptic-function potential. Finally, in 
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Sec. IV we make some concluding remarks and describe the future objectives of our research 



program. 

II. PREVIOUS RESULTS ON COMPLEX CLASSICAL MECHANICS 

During the past decade there has been an active research program to extend quantum 
mechanics into the complex domain. Specifically, it has been shown that the requirement 
that a Hamiltonian be Dirac Hermitian (we say that a Hamiltonian is Dirac-Hermitian if 
H = H\ where f represents the combined operations of complex conjugation and matrix 
transposition) may be broadened to include complex non-Dirac-Hermitian Hamiltonians that 
are VT symmetric. This much wider class of Hamiltonians are physically acceptable because 
they possess two crucial features: (i) their eigenvalues are all real, and (ii) they describe 
unitary time evolution. (We say that a Hamiltonian is VT symmetric if it is invariant under 
combined spatial reflection and time reversal.) 

An example of a class of Hamiltonians that is not Dirac Hermitian but which is VT 
symmetric is given by 

H = p'^ + x\ixY (e>0). (1) 

When e = O, H, which represents the familiar quantum harmonic oscillator, is Dirac Hermi- 
tian. While H is no longer Dirac Hermitian when e increases from O, H continues to be VT 
symmetric, and its eigenvalues continue to be real, positive, and discrete P-l9]. Because a 
PT-symmetric quantum system in the complex domain retains the fundamental properties 
required of a physical quantum theory, much theoretical research on such systems has been 
published and recent experimental observations have confirmed some theoretical predictions 
[10HI3]. 

Complex quantum mechanics has proved to be so interesting that the research activity on 
VT quantum mechanics has motivated studies of complex classical mechanics. In the study 
of complex systems the complex as well as the real solutions to Hamilton's differential equa- 
tions of motion are considered. In this generalization of conventional classical mechanics, 
classical particles are not constrained to move along the real axis and may travel through 
the complex plane. 

Early work on the particle trajectories in complex classical mechanics is reported in 
Refs. [21 [T^. Subsequently, detailed studies of the complex extensions of various one- 
dimensional conventional classical-mechanical systems were undertaken: The remarkable 
properties of complex classical trajectories are examined in Refs. [T3] - [T^ . Higher dimen- 
sional complex classical-mechanical systems, such as the Lotka-Volterra equations for popu- 
lation dynamics and the Euler equations for rigid body rotation are discussed in Refs. pU] . 
The complex PT-symmetric Korteweg-de Vries equation has also been studied [2T1 - I27] . 

The objective in extending classical mechanics into the complex domain is to enhance 
our understanding of subtle mathematical and physical phenomena. For example, it was 
found that some of the complicated properties of chaotic systems become more transparent 
when extended into the complex domain |28] . Also, studies of exceptional points of complex 
systems have revealed interesting and potentially observable effects [29l 130] . Finally, recent 
work on the complex extension of quantum probability density constitutes an advance in 
our understanding of the quantum correspondence principle pi] . 

An elementary example that illustrates the extension of a conventional classical- 
mechanical system into the complex plane is given by the classical harmonic oscillator. 
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FIG. 1: Classical trajectories in the complex plane for the harmonic-oscillator Hamiltonian 
H = + x'^. These trajectories are nested ellipses. Observe that when the harmonic oscilla- 
tor is extended into the complex-x domain, the classical particles may pass through the classically 
forbidden regions on the real axis outside the turning points. When the trajectories cross the real 
axis, they are orthogonal to it. 



whose Hamiltonian is given in ([T]) with e = 0. The Standard classical equations of motion 
for this system are 

X = 2p, p = —2x. (2) 

However, we now treat the coordinate variable z{t) and the momentum variable p{t) to be 
complex functions of time t. That is, we consider this system to have one complex degree of 
freedom. Thus, the equations of motion become 

r = 2u, s = 2v, ii = —2r, v = —2s, (3) 

where the complex coordinate is x = r + is and the complex momentum is p = u + iv. For 
a particle having real energy E and initial position r(0) = a > a/E, s(0) = O, the solution 
to Q is 

r(t) = acos(2t), s{t) = Va^ -Esin{2t). (4) 

Thus, the possible classical trajectories are a family of ellipses parametrized by the initial 
position a: 

Five of these trajectories are shown in Fig. [Tj Each trajectory has the same period T = n. 
The degenerate ellipse, whose foci are the turning points at x = ±\/E, is the familiar real 
solution. Note that classical particles may visit the real axis in the classically forbidden 
regions |x| > VE, but that the elliptical trajectories are orthogonal rather than parallel to 
the real-x axis. 

In general, when classical mechanics is extended into the complex domain, classical par- 
ticles are allowed to enter the classically forbidden region. However, in the forbidden region 
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there is no particle flow parallel to the real axis and the flow of classical particles is orthogo- 
nal to the axis. This feature is analogous to the vanishing flux of energy in the case of total 
internal refiection. 

In the case of total internal refiection when the angle of incidence is less than a critical 
value, there is a refiected wave but no transmitted wave. The electromagnetic field does 
cross the boundary and this field is attenuated exponentially in a few wavelengths beyond 
the interface. Although the field does not vanish in the classically forbidden region, there 
is no fiux of energy; that is, the Poynting vector vanishes in the classically forbidden region 
beyond the interface. We emphasize that in the physical world the cutoff at the boundary 
between the classically allowed and the classically forbidden regions is not perfectly sharp. 
For example, in classical optics it is known that below the surface of an imperfect conductor, 
the electromagnetic fields do not vanish abruptly. Rather, they decay exponentially as 
functions of the penetration depth. This effect is known as skin depth [52] . 

Another model that illustrates the properties of the classically allowed and classically 
forbidden regions is the anharmonic oscillator, whose Hamiltonian is 

H = \p^ + (6) 

For this Hamiltonian there are four turning points, two on the real axis and two on the 
imaginary axis. When the energy is real and positive, all the classical trajectories are closed 
and periodic except for two special trajectories that begin at the turning points on the 
imaginary axis. Four trajectories for the case E = 1 are shown in Fig. |2} 

The topology of the classical trajectories changes dramatically if the classical energy is 
allowed to be complex: When ImE' 7^ O, the classical paths no longer closed. This feature 
is illustrated in Fig. 3, which shows the path of a particle of energy £^ = 1 + O.li in an 
anharmonic potential. 

The observation that classical orbits are closed and periodic when the energy is real and 
open and nonperiodic when the energy is complex was made in Ref. [20] and studied in detail 
in Ref. [33j- In these references it is emphasized that the Bohr-Sommerfeld quantization 
condition 

dxp = [n+ ^) n (7) 

can only be applied if the classical orbits are closed. Thus, there is a deep connection 
between real classical energies and the existence of associated real quantum eigenvalues. 

It was further argued in Ref. ^3] that the measurement of a quantum energy is inherently 
imprecise because of the time-energy uncertainty principle AE At > h/2. Specifically, since 
there is not an infinite amount of time in which to make a quantum energy measurement, 
we expect that the uncertainty in the energy AE is nonzero. If we then suppose that 
this uncertainty has an imaginary component, it follows that in the corresponding classical 
theory, while the particle trajectories are almost periodic, the orbits do not close exactly. 
The fact that the classical orbits with complex energy are not closed means that in complex 
classical mechanics one can observe tunneling-like phenomena that one normally expects to 
find only in quantum systems. 

We illustrate such tunneling-like phenomena by considering a particle in a quartic double- 
well potential V{x) = x^ — 5x'^. Figure |4] shows eight possible complex classical trajectories 
for a particle of real energy E = —1. Each of these trajectories is closed and periodic. 
Observe that for this energy the trajectories are localized either in the left well or the right 
well and that no trajectory crosses from one side to the other side of the imaginary axis. 
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FIG. 2: Classical trajectories x{t) in the complex-x plane for the anharmonic-oscillator Hamiltonian 
H = ^p^ + Ali trajectories represent a particle of energy E = 1. There is one real trajectory 
that oscillates between the turning points at x = ±1 and an infinite family of nested complex 
trajectories that enclose the real turning points but lie inside the imaginary turning points at 
ibi. (The turning points are indicated by dots.) Two other trajectories begin at the imaginary 
turning points and drift off to infinity along the imaginary- 3; axis. Apart from the trajectories 
beginning at ±i, all trajectories are closed and periodic. Ali orbits in this figure have the same 
period v^vr/^r (|) /T (|) = 3.70815 . . .. 

What happens if we allow the classical energy to be complex? In this case the classical 
trajectory is no longer closed. However, it does not spiral out to infinity like the trajectory 
shown in Fig. 3. Rather, the trajectory in Fig. [5] unwinds around a pair of turning points 
for a characteristic length of time and then crosses the imaginary axis. At this point the 
trajectory does something remarkable: Rather than continuing its outward journey, it spirals 
inward towards the other pair of turning points. Then, never crossing itself, the trajectory 
turns outward again, and after the same characteristic length of time, returns to the vicinity 
of the first pair of turning points. This oscillatory behavior, which shares the qualitative 
characteristics of strange attractors, continues forever but the trajectory never crosses itself. 
As in the case of quantum tunneling, the particle spends a long time in proximity to a 
given pair of turning points before crossing the imaginary axis to the other pair of turning 
points. On average, the classical particle spends equal amounts of time on either side of 
the imaginary axis. Interestingly, we find that as the imaginary part of the classical energy 
increases, the characteristic "tunneling" time decreases in inverse proportion, just as one 
would expect of a quantum particle. 

Having described the tunneling-like behavior of a classical particle having complex energy 
in a double well, we examine the case of such a particle in a periodic potential. Physically, 
this corresponds to a classical particle in a crystal lattice. A simple physical system that 
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FIG. 3: A single classical trajectory in the complex-j; plane for a particle governed by the 
anharmonic-oscillator Hamiltonian H = + x^. This trajectory begins at x = 1 and repre- 
sents the complex path of a particle whose energy E = 1 + O.li is complex. The trajectory is not 
closed or periodic. The four turning points are indicated by dots. The trajectory does not cross 
itself. 
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FIG. 4: Eight classical trajectories in the complex-a; plane representing a particle of energy E = — 1 
in the potential x^ — f)x'^. The turning points are located at x = ±2.19 and x = ±0.46 and are 
indicated by dots. Because the energy is real, the trajectories are all closed. The classical particle 
stays in either the right-half or the left-half plane and cannot cross the imaginary axis. Thus, when 
the energy is real, there is no effect analogous to tunneling. 

has a periodic potential consists of a simple pendulum in a uniform gravitational field |31] ■ 
Consider a pendulum consisting of a bob of mass m and a string of length L in a uniform 
gravitational field of magnitude g (see Fig. |6|. The gravitational potential energy of the 
system is defined to be zero at the height of the pivot point of the string. The pendulum 
bob swings through an angle 9. Therefore, the horizontal and vertical cartesian coordinates 



7 




Re X 



FIG. 5: Classical trajectory of a particle moving in the complex-x plane under the influence of a 
double-well x^ — potential. The particle has complex energy E = —l — i and thus its trajectory 
does not close. The trajectory spirals outward around one pair of turning points, crosses the 
imaginary axis, and then spirals inward around the other pair of turning points. It then spirals 
outward again, crosses the imaginary axis, and goes back to the original pair of turning points. 
The particle repeats this behavior endlessly but at no point does the trajectory cross itself. This 
classical-particle motion is analogous to the behavior of a quantum particle that repeatedly tunnels 
between two classically allowed regions. Here, the particle does not disappear into the classically 
forbidden region during the tunneling process; rather, it moves along a well-defined path in the 
complex-x plane from one well to the other. 



X and Y are X = Lsind and Y = —LcosO, which gives velocities X = LOsmO and 
Y = —L9cos9. The potential and kinetic energies are V = —mgLcosO and T = ^m(X^ + 
Y"^) = ^mL^9'^. The Hamiltonian H = T + V for the pendulum is therefore 

H = \mL'^e^ -mgLcosO. (8) 

Without loss of generality we set m = 1, (? = 1, and L = 1 and then make the change of 
variable 9 ^ x to get 

H = — cosx, (9) 
where p = x. The classical equations of motion for this Hamiltonian are 

dH dH ,^ , 

X = — — = p, p = — — — = — smx. (10) 
op ox 
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FIG. 6: Configuration of a simple pendulum of mass m in a uniform gravitational field of strength 
g. The length of the string is L. The pendulum swings through an angle 6. We define the potential 
energy to be O at the height of the pivot. 
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FIG. 7: Classical trajectories in the complex-x plane for a particle of energy E = —0.09754 in a 
— cosx potential. The motion is periodic and the particle remains confined to a cell of width 2tt. 
Five trajectories are shown for each cell. The trajectories shown here are the periodic analogs of 
the trajectories shown in Figs. [T] and [2| 

The Hamiltonian H for this system is a constant of the motion and thus the energy E' is a 
time-independent quantity. 

If we take the energy E to be real, we find that the classical trajectories are confined to 
cells of horizontal width 27r, as shown in Fig. [7j This is the periodic analog of Figs. [T] and |2} 

If the energy of the classical particle in a periodic potential is taken to be complex, the 
particle begins to hop from well to well in analogy to the behavior of the particle in Fig. |5| 
This hopping behavior is displayed in Fig. [8| 

The most interesting analogy between quantum mechanics and complex classical me- 
chanics is established by showing that there exist narrow conduction bands in the periodic 
potential for which the quantum particle exhibits resonant tunneling and the complex clas- 




FIG. 8: A tunneling trajectory for the Hamiltonian ^ with E = 0.1 — 0.1 5i. The classical particle 
hops from well to well in a random-walk fashion. The particle starts at the origin and then hops 
left, right, left, left, right, left, left, right, right. This is the sort of behavior normally associated 
with a particle in a crystal at an energy that is not in a conduction band. At the end of this 
simulation the particle is situated to the left of its initial position. The trajectory never crosses 
itself. 



sical particle exhibits unidirectional hopping [33|l35]. This qualitative behavior is illustrated 
in Fig. 1} 

A detailed numerical analysis shows that the classical conduction bands have a narrow 



but finite width (see Fig. 10). Two magnified portions of the conduction bands in Fig. 10 



are shown in Fig. VI. These magnifications show that the edges of the conduction bands are 
sharply defined. 



III. CLASSICAL PARTICLE IN A COMPLEX ELLIPTIC POTENTIAL 

Having reviewed in Sec. |Tl]the behavior of complex classical trajectories for trigonometric 
potentials, in this section we give a brief glimpse of the rich and interesting behavior of clas- 
sical particles moving in elliptic potentials. Elliptic potentials are natural doubly-periodic 
generalizations of trigonometric potentials. 

The Hamiltonian that we have chosen to study is a simple extension of that in ([9]): 

H = -Cn{x,k), (11) 

where Cn{x,k) is a cnoidal function [361 EZ]- When the parameter = O, the cnoidal 
function reduces to the singly periodic function cosx and when k = 1 the cnoidal function 
becomes tanhx. When O < A; < 1, the cnoidal function is periodic in both the real and 
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FIG. 9: A classical particle exhibiting a behavior analogous to that of a quantum particle in a 
conduction band that is undergoing resonant tunneling. Unlike the particle in Fig. [8j this classical 
classical particle tunnels in one direction only and drifts at a constant average velocity through the 
potential. 



imaginary directions and it is meromorphic (analytic in the finite-x plane except for pole 
singularities) and has infinitely many double poles. The real part of the cnoidal potential 



Cn(a;, k) is plotted in Fig. 12 



The classical particle trajectories satisfy Hamilton's equations 



X 



P 



dH _ 
dp 
dH 

dx 



P, 



Sn{x, fc)Dn(x, k). 



(12) 



The trajectories for k > O are remarkable in that the classical particles seem to prefer to 



move vertically rather than horizontally. In Fig. 13 a trajectory, similar to that in Fig. [8], 
is shown for the case k = 0. This trajectory is superimposed on a plot of the real part of 
the cosine potential. The particle oscillates horizontally. In Fig. [14] a complex trajectory for 
the case k = 1/5 is shown. This more exotic path escapes from the initial pair of turning 
points, and rather than "tunneling" to a horizontally adjacent pair of turning points, it 
travels downward. The ensuing wavy vertical motion passes close to many poles before 
the particle gets captured by another pair of turning points. The particle winds inwards 
and outwards around these turning points and eventually returns to the original pair of 
turning points. After escaping from these turning points again, the particle now moves in 
the positive-imaginary direction. 



-1.0 -0.5 0.0 0.5 1.0 

Re[E] 

FIG. 10: Complex-energy plane showing those energies that lead to tunneling (hopping) behavior 
and those energies that give rise to conduction. Hopping behavior is indicated by a hyphen - and 
conduction is indicated by an X. The symbol & indicates that no tunnehng takes place; tunnehng 
does not occur for energies whose imaginary part is close to 0. In some regions of the energy plane 
we have done very intensive studies and the X's and -'s are densely packed. This picture suggests 
the features of band theory: If the imaginary part of the energy is taken to be —0.9, then as the real 
part of the energy increases from —1 to +1, five narrow conduction bands are encountered. These 
bands are located near Re i? = —0.95, —0.7, —0.25, 0.15, 0.7. This picture is symmetric about 
luiE = O and the bands get thicker as llmi?! increases. A total of 68689 points were classified 
to make this plot. In most places the resolution (distance between points) is 0.01, but in several 
regions the distance between points is shortened to 0.001. The regions indicated by arrows are 



blown up in Fig. 11 



It is clear that classical trajectories associated with doubly periodic potentials have an im- 
mensely interesting structure and should be investigated in much greater detail to determine 



if there is a behavior analogous to band structure shown in Figs. 10 and 11 



IV. SUMMARY AND DISCUSSION 



The relationship between quantum mechanics and classical mechanics is subtle. Quantum 
mechanics is essentially wavelike; probability amplitudes are described by a wave equation 
and physical observations involve such wavelike phenomena as interference patterns and 
nodes. In contrast, classical mechanics describes the motion of particles and exhibits none of 
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FIG. 11: Detailed portions of the complex-energy plane shown in Fig. [TO] containing a conduction 
band. Note that the edge of the conduction band, where tunnehng (hopping) behavior changes 
over to conducting behavior, is very sharp. 




FIG. 12: Real part of the cnoidal elhptic function Cn(x,A;) in the complex-rc plane for two values 
of k, k = O and k = 1/5. When k = O, this cnoidal function reduces to the trigonometric function 
cosx, and this function grows exponentially in the imaginary-x direction. When A: > O, the cnoidal 
function is doubly periodic; that is, periodic in the real- a; and in the imaginary-x directions. While 
cosx is entire, the cnoidal functions for 7^ O are meromorphic and have periodic double poles. 

these wavelike features. Nevertheless, there is a deep connection between quantum mechanics 
and complex classical mechanics. In the complex domain the classical trajectories exhibit a 
remarkable behavior that is analogous to quantum tunneling. 

Periodic potentials exhibit a surprising and intricate feature that closely resembles quan- 
tum band structure. It is especially noteworthy that the classical bands, just like the quan- 
tum bands, have finite width. 

Our early work on singly periodic potentials strongly suggests that further detailed analy- 
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FIG. 13: Complex classical trajectory for a particle of energy E = 0.5 + 0.05i in a cosine potential 
(a cnoidal potential with /c = 0). The trajectory begins at the red dot at x = 0.5i on the left side of 
the figure and spirals outward around the left pair of turning points. The trajectory then "tunnels" 
to the right and spirals inward and then outward around the right pair of turning points. In the 



background is a plot of the real part of the cosine potential, which is shown in detail in Fig. 12 



sis should be done on doubly periodic potentials. Doubly periodic potentials are particularly 
interesting because they have singularities. Two important and so far unanswered questions 
are as foUows: (i) Does a complex classical particle in a doubly periodic potential undergo 
a random walk in two dimensions and eventually visit all lattice sites? (ii) Are there spe- 
cial bands of energy for which the classical particle no longer undergoes random hopping 
behavior and begins to drift in one direction through the lattice? 
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FIG. 14: Complex classical trajectory for a particle of energy E = 0.5 + 0.05z in a cnoidal potential 
with k = 1/5. The trajectory begins at the red dot at a; = 0.5i near the top of the central pane. 
The trajectory spirals outward around the two turning pomts and, when it gets very close to a 
pole, it suddenly begins to travel downward in a wavy fashion in the negative-imaginary direction. 
The trajectory bobs and weaves past many double poles and continues on in the left pane of the 
figure. It eventually gets very close to a double pole, gets trapped, and spirals inward towards a 
pair of turning points. After spiraling inwards, it then spirals outwards (never crossing itself) and 
goes upward along a path extremely close to the downward wavy path. It is then recaptured by 
the original pair of turning points in the central pane. After spiraling inwards and outwards once 
more it now escapes and travels upward, bobbing and weaving along a wavy path in the right pane 
of the figure. Evidently, the particle trajectory strongly prefers to move vertically upward and 
downward, and not horizontally. The vertical motion distinguishes cnoidal trajectories from those 



in Figs. ^ M and 13 associated with the cosine potential. 



was used to generate the figures in this paper. 
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